Genetic (co-)variation of methane emissions, efficiency, and production traits in Danish Holstein cattle along and across lactations.

Methane emissions in ruminant livestock has become a hot topic, given the pressure to reduce greenhouse gas emissions drastically in the European Union over the next 10 to 30 yr. During the 2021 United Nations Climate Change conference, countries also made collective commitments to curb methane emissions by 2050. Genetic selection for low-methane-emitting animals, particularly dairy cows, is one possible strategy for mitigation. However, it is essential to understand how methane emissions in lactating animals vary along lactation and across lactations. This understanding is useful when making decisions for future phenotyping strategies, such as the frequency and duration of phenotyping within and across lactations. Therefore, the objectives of this study were to estimate (1) genetic parameters for 2 methane traits: methane concentration (MeC) and methane production (MeP) at 2 parity levels in Danish Holstein cows across the entire lactation using random regression models; (2) genetic correlations within and between methane traits across the entire lactation; and (3) genetic correlations between the methane traits and economically important traits throughout first lactation. Methane concentration (n = 19,639) records of 575 Danish Holstein cows from a research farm measured between 2013 and 2020 were available. Subsequently, CH4 production in grams/day (MeP; n = 13,866) was calculated; MeP and MeC for first and second lactation (L1 and L2) were analyzed as separate traits: MeC_L1, MeP_L1, MeC_L2, and MeP_L2. Heritabilities, variance components, and genetic correlations within and between the 4 CH4 traits were estimated using random regression models with Legendre polynomials. The additive genetic and permanent environmental effects were modeled using second-order Legendre polynomial for lactation weeks. Estimated heritabilities for MeP_L1 ranged between 0.11 and 0.49, for MeC_L1 between 0.10 and 0.28, for MeP_L2 between 0.14 and 0.36, and for MeC_L2 between 0.13 and 0.29. In general, heritability estimates of MeC traits were lower and more stable throughout lactation and were similar between lactations compared with MeP. Genetic correlations (within trait) at different lactation weeks were generally highly positive (0.7) for most of the first lactation, except for the correlation of early lactation (<10 wk) with late lactation (>40 wk) where the correlation was the lowest (<0.5). Genetic correlations between methane traits were moderate to highly correlated during early and mid lactation. Finally, MeP_L1 has stronger genetic correlations with energy-corrected milk and dry matter intake compared with MeC_L1. In conclusion, both traits are different along (and across) lactation(s) and they correlated differently with production, maintenance, and intake traits, which is important to consider when including one of them in a future breeding objective.

Variation in volatile organic compounds in Atlantic salmon mucus is associated with resistance to salmon lice infection.

Salmon lice are ectoparasites that threaten wild and farmed salmonids. Artificial selection of salmon for resistance to the infectious copepodid lice stage currently relies on in vivo challenge trials on thousands of salmon a year. We challenged 5750 salmon with salmon lice (Lepeophtheirus salmonis) from two distinct farmed strains of salmon in two separate trials. We found that volatile organic compounds (VOC), 1-penten-3-ol, 1-octen-3-ol and 6-methyl-5-hepten-2-one in the mucus of the salmon host after salmon lice infection, were significantly associated with lice infection numbers across a range of water temperatures (5 °C, 10 °C, 17 °C). Some VOCs (benzene, 1-octen-3-ol and 3,5,5-trimethyl-2-hexene) were significantly different between lines divergently selected for salmon lice resistance. In a combined population assessment, selected VOCs varied between families in the range of 47- 59% indicating a genetic component and were positively correlated to the salmon hosts estimated breeding values 0.59-0.74. Mucosal VOC phenotypes could supplement current breeding practices and have the potential to be a more direct and ethical proxy for salmon lice resistance provided they can be measured prior to lice infestation.

Genome-wide association study for methane emission traits in Danish Holstein cattle.

Selecting for lower methane emitting cows requires insight into the most biologically relevant phenotypes for methane emission, which are close to the breeding goal. Several methane phenotypes have been suggested over the last decade. However, the (dis)similarity of their underlying genetic architecture and correlation structures are poorly understood. Therefore, the objective of this study was to test association of SNP and genomic regions through GWAS on 8 CH4 emission traits in Danish Holstein cattle. The traits studied were methane concentration (MeC; ppm), methane production (MeP ; g/d), 2 definitions of residual methane (RMETc and RMETp: MeC and MeP regressed on metabolic body weight and energy-corrected milk, respectively), 2 definitions of methane intensity (MeI; MeIc = MeC/ECM and MeIp = MeP/ECM); 2 definitions of methane yield per kilogram of dry matter intake (MeY; MeYc = MeC/dry matter intake and MeYp = MeP/dry matter intake). A total of 1,962 cows with genotypes (Illumina BovineSNP50 Chip or Eurogenomic custom SNP chip) and repeated records of the above-mentioned 8 methane traits were analyzed. Strong associations were found with 3 traits (MeC, MeP, and MeYc) on chromosome 13 and with 5 traits (MeC, MeP, MeIp, MeYp, and MeYc) on chromosome 26. For MeIc, MeIp, RMETc, MeYc, and MeYp, some suggestive association signals were identified on chromosome 1. Genomic segments of 1 Mbp (n = 2,525) were tested for their association with these traits, which identified between 33 to 54 significantly associated regions. In a pairwise comparison, MeC and MeP were the traits that shared the highest number of significant segments (17). The same trend was observed when comparing SNP significantly associated with the traits MeC and MeP shared from 23 to 25 SNP (most of which were located in chromosomes 11, 13, and 26). Based on our results on GWAS and genetic correlations, we conclude that MeC is (genetically) more closely linked to MeP than any of the other methane traits analyzed.

Breeding for reduced methane emission and feed-efficient Holstein cows: An international response.

Selecting for lower methane (CH4) emitting animals is one of the best approaches to reduce CH4 given that genetic progress is permanent and cumulative over generations. As genetic selection requires a large number of animals with records and few countries actively record CH4, combining data from different countries could help to expedite accurate genetic parameters for CH4 traits and build a future genomic reference population. Additionally, if we want to include CH4 in the breeding goal, it is important to know the genetic correlations of CH4 traits with other economically important traits. Therefore, the aim of this study was first to estimate genetic parameters of 7 suggested methane traits, as well as genetic correlations between methane traits and production, maintenance, and efficiency traits using a multicountry database. The second aim was to estimate genetic correlations within parities and stages of lactation for CH4. The third aim was to evaluate the expected response of economically important traits by including CH4 traits in the breeding goal. A total of 15,320 methane production (MeP, g/d) records from 2,990 cows belonging to 4 countries (Canada, Australia, Switzerland, and Denmark) were analyzed. Records on dry matter intake (DMI), body weight (BW), body condition score, and milk yield (MY) were also available. Additional traits such as methane yield (MeY; g/kg DMI), methane intensity (MeI; g/kg energy-corrected milk), a genetic standardized methane production, and 3 definitions of residual methane production (g/d), residual feed intake, metabolic BW (MBW), BW change, and energy-corrected milk were calculated. The estimated heritability of MeP was 0.21, whereas heritability estimates for MeY and MeI were 0.30 and 0.38, and for the residual methane traits heritability ranged from 0.13 to 0.16. Genetic correlations between different methane traits were moderate to high (0.41 to 0.97). Genetic correlations between MeP and economically important traits ranged from 0.29 (MY) to 0.65 (BW and MBW), being 0.41 for DMI. Selection index calculations showed that residual methane had the most potential for inclusion in the breeding goal when compared with MeP, MeY, and MeI, as residual methane allows for selection of low methane emitting animals without compromising other economically important traits. Inclusion of residual feed intake in the breeding goal could further reduce methane, as the correlation with residual methane is moderate and elicits a favorable correlated response. Adding a negative economic value for methane could facilitate a substantial reduction in methane emissions while maintaining an increase in milk production.

Multitrait genomic prediction of methane emissions in Danish Holstein cattle.

In dairy cattle, selecting for lower methane-emitting animals is one of the new challenges of this decade. However, genetic selection requires a large number of animals with records to get accurate estimated breeding values (EBV). Given that CH4 records are scarce, the use of information on routinely recorded and highly correlated traits with CH4 has been suggested to increase the accuracy of genomic EBV (GEBV) through multitrait (genomic) prediction. Therefore, the objective of this study was to evaluate accuracies of prediction of GEBV for CH4 by including or omitting CH4, energy-corrected milk (ECM), and body weight (BW) as well as genotypic information in multitrait analyses across 2 methods: BLUP and single-step genomic BLUP (SSGBLUP). A total of 2,725 cows with CH4 concentration in breath (14,125 records), BW (61,667 records), and ECM (61,610 records) were included in the analyses. Approximately 2,000 of these cows were genotyped or imputed to 50K. Ten cross-validation groups were formed by randomly grouping paternal half-sibs. Five scenarios were performed: (1) base scenario with only CH4 information; (2) without CH4, but with information from BW, ECM, or BW+ECM only in reference population; (3) without CH4, but with information from BW, ECM, or BW+ECM in both validation and reference population; (4) with CH4 information and BW, ECM, or BW+ECM information only in the reference population; and (5) with CH4 information and BW, ECM, or BW+ECM information in both validation and reference population. As a result, for each method (BLUP, SSGBLUP), 13 sub-scenarios were performed, 1 from scenario 1, and 3 for each of the subsequent 4 scenarios. The average accuracy of GEBV for CH4 in the base scenario was 0.32 for BLUP and 0.42 for SSGBLUP, and it ranged from 0.10 in scenario 2 to 0.78 in scenario 5 across methods. In terms of bias, the base scenario 1 was unbiased for SSGBLUP; similar results were achieved with scenario 5. Including information on ECM increased the accuracy of GEBV for CH4 by up to 61%, whereas adding information on both traits (BW and ECM) increased the accuracy by up to 90%. Scenarios that did not include CH4 in the reference population had the lowest correlations (0.17-0.33) with single-trait CH4 GEBV, and scenarios with CH4 in the reference population had the highest correlations (0.41-0.81). Thus, failure to include CH4 in future reference populations results in predicted CH4 GEBV, which cannot be used in practical selection. Therefore, recording CH4 in more animals remains a priority. Finally, multiple-trait genomic prediction using routinely recorded BW and ECM leads to higher prediction accuracies than traditional single-trait genomic prediction for CH4 and is a viable solution for increasing the accuracies of GEBV for scarcely recorded CH4 in practice.

Review: Genetic and genomic selection as a methane mitigation strategy in dairy cattle.

Over the last decade, extensive research effort has been placed on developing methane mitigation strategies in ruminants. Many disciplines on animal science disciplines have been involved, including nutrition and physiology, microbiology and genetic selection. To date, few of the suggested strategies have been implemented because: (1) methane emissions currently have no direct or indirect economic value for farmers, with no financial incentive to change practices and (2) most strategies have limited, or no, long-term effects. Consequently, there is a fundamental need for research on methane mitigation strategies across disciplines. Coordinated international initiatives similar to METHAGENE could represent highly relevant coordination tool of collaboration between countries, facilitating knowledge exchange, sharing concerns and building future collaborations.

Phenotypic modeling of residual feed intake using physical activity and methane production as energy sinks.

Residual feed intake (RFI) is a measure of feed efficiency in dairy cattle. This study modeled phenotypic RFI of first- and second-parity Holstein and Jersey dairy cows within 9 lactation segments (consecutive segments of 4 wk each) covering the first 36 lactation weeks. We aimed to evaluate physical activity and daily methane production as additional energy sinks in the estimation of RFI, to examine the correlations of RFI among the first 36 wk of lactation (WOL), and to evaluate whether parities and breeds show similar results. Records for first-parity Holstein (n = 449), second-parity Holstein (n = 298), first-parity Jersey (n = 195), and second-parity Jersey cows (n = 146) were used. Model 1 included the following energy sinks: energy-corrected milk yield, metabolic body weight (BW), body condition score (BCS), daily changes in BW (ΔBW) and BCS (ΔBCS), and physical activity. Model 2 was based on a subset of the data and only for Holstein cows, and included the same energy sinks as Model 1, plus daily methane production. The trajectories of segment-specific partial regression coefficients (PRC) of DMI on activity were similar across parities but differed slightly between breeds. For daily methane production, the trajectory in PRC decreased over lactation segments for first- and second-parity Holstein cows. The trajectories in PRC of DMI on energy-corrected milk yield, metabolic BW, BCS, and ΔBW were generally similar across parities, except for ΔBCS. Activity accounted for on average 7.3, 6.8, 7.2, and 6.4% of DMI for first-parity Holsteins, second-parity Holsteins, first-parity Jerseys, and second-parity Jerseys, respectively. Methane losses accounted for 8.7% and 8.5% of DMI for first- and second-parity Holstein cows, respectively. Repeatability estimates for RFI over 36 WOL for Model 1 were 0.63 for first-parity Holsteins, 0.65 for second-parity Holsteins, 0.76 for first-parity Jerseys, and 0.80 for second-parity Jerseys. For Model 2, the estimates were 0.59 and 0.61 for first- and second-parity Holstein cows, respectively. Correlations of RFI between WOL varied in strength, with weak correlations for the first 2 to 3 WOL with other WOL. In conclusion, physical activity and daily methane production accounted for part of DMI, and RFI of dairy cattle is not the same trait throughout lactation.

Can greenhouse gases in breath be used to genetically improve feed efficiency of dairy cows?

There is considerable interest in improving feed utilization of dairy cattle while limiting losses to the environment (i.e., greenhouse gases, GHG). To breed for feed-efficient or climate-friendly cattle, it is first necessary to obtain accurate estimates of genetic parameters and correlations of feed intake, greenhouse gases, and production traits. Reducing dry matter take (DMI) requirements while maintaining production has high economic value to farmers, but DMI is costly to record and thus limited to small research or nucleus herds. Conversely, enteric methane (CH4) currently has no economic value, is also costly to record, and is limited to small experimental trials. However, breath gas concentrations of methane (CH4c) and carbon dioxide (CO2c) are relatively cheap to measure at high throughput under commercial conditions by installing sniffers in automated milking stations. The objective of this study was to assess the genetic correlations between DMI, body weight (BW), fat- and protein-corrected milk yield (FPCM), and GHG-related traits: CH4c and CO2c from Denmark (DNK) and the Netherlands (NLD). A second objective was to assess the genetic potential for improving feed efficiency and the added benefits of using CH4c and CO2c as indicators. Feed intake data were available on 703 primiparous cows in DNK and 524 in NLD; CH4c and CO2c records were available on 434 primiparous cows in DNK and 656 in NLD. The GHG-related traits were heritable (e.g., CH4c h2: DNK = 0.26, NLD = 0.15) but were differentially genetically correlated with DMI and feed efficiency in both magnitude and sign, depending on the population and the definition of feed efficiency. Across feed efficiency traits and DMI, having bulls with 100 daughters with FPCM, BW, and GHG traits resulted in sufficiently high accuracy to almost negate the need for DMI records. Despite differences in genetic correlation structure, the relatively cheap GHG-related traits showed considerable potential for improving the accuracy of breeding values of highly valuable feed intake and feed efficiency traits.

Is rumination time an indicator of methane production in dairy cows?

As long as large-scale recording of expensive-to-measure and labor-consuming traits, such as dry matter intake (DMI) and CH4 production (CH4P), continues to be challenging in practical conditions, alternative traits that are already routinely recorded in dairy herds should be investigated. An ideal indicator trait must, in addition to expressing genetic variation, have a strong correlation with the trait of interest. Our aim was to estimate individual level and phenotypic correlations between rumination time (RT), CH4P, and DMI to determine if RT could be used as an indicator trait for CH4P and DMI. Data from 343 Danish Holstein cows were collected at the Danish Cattle Research Centre for a period of approximately 3 yr. The data set consisted of 14,890 records for DMI, 15,835 for RT, and 6,693 for CH4P. Data were divided in primiparous cows only (PC) and all cows (MC), and then divided in lactation stage (early, mid, late, and whole lactation) to analyze the changes over lactation. Linear mixed models, including an animal effect but no pedigree, were used to estimate the correlations among traits. Phenotypic and individual level correlations between RT and both CH4P and DMI were close to zero, regardless of lactation stage and data set (PC or MC). However, CH4P and DMI were highly correlated, both across lactation stages and data sets. In conclusion, RT is unsuitable to be used as an indicator trait for either CH4P or DMI. Our study failed to validate RT as a useful indicator trait for both CH4P and DMI, but more studies with novel phenotypes can offer different approaches to select and incorporate important yet difficult to record traits into breeding goals and selection indexes.

Review: Selecting for improved feed efficiency and reduced methane emissions in dairy cattle.

It may be possible for dairy farms to improve profitability and reduce environmental impacts by selecting for higher feed efficiency and lower methane (CH4) emission traits. It remains to be clarified how CH4 emission and feed efficiency traits are related to each other, which will require direct and accurate measurements of both of these traits in large numbers of animals under the conditions in which they are expected to perform. The ranking of animals for feed efficiency and CH4 emission traits can differ depending upon the type and duration of measurement used, the trait definitions and calculations used, the period in lactation examined and the production system, as well as interactions among these factors. Because the correlation values obtained between feed efficiency and CH4 emission data are likely to be biased when either or both are expressed as ratios, therefore researchers would be well advised to maintain weighted components of the ratios in the selection index. Nutrition studies indicate that selecting low emitting animals may result in reduced efficiency of cell wall digestion, that is NDF, a key ruminant characteristic in human food production. Moreover, many interacting biological factors that are not measured directly, including digestion rate, passage rate, the rumen microbiome and rumen fermentation, may influence feed efficiency and CH4 emission. Elucidating these mechanisms may improve dairy farmers ability to select for feed efficiency and reduced CH4 emission.